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In  phosphoric  acid  fuel  cell  (PAFC)  stack,  reduction  of  the  phosphoric  acid  that  is  impregnated  in  the  cell  is 
a  major  factor  to  restrict  the  operating  life.  In  this  paper,  the  phosphoric  acid  reduction  by  the  evaporative 
and  condensational  dissipation  was  evaluated  by  numerical  analysis.  The  calculations  that  include  the 
behavior  of  heat  transfer,  gas  flow,  and  evaporation  and  condensation  of  phosphoric  acid  were  conducted 
for  the  model  cell  with  conditions  that  correspond  to  an  on  site  stack.  The  phosphoric  acid  was  considered 
as  a  composite  of  phosphorus  pentoxide  and  water,  and  the  evaporation  and  condensation  rates  were 
derived  based  on  the  nucleation  theory  for  a  two  components  system.  The  phosphoric  acid  distributions 
in  the  vapor  phase  and  liquid  phase  at  the  process  gas  and  electrode  were  calculated  for  the  duration 
of  3300  h,  and  as  a  result  of  those,  the  exhausting  and  remaining  amount  of  phosphoric  acid  for  the  cell 
were  evaluated.  The  analysis  results  were  compared  with  the  experimental  results  for  the  model  cell. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  phosphoric  acid  fuel  cell  (PAFC)  stack,  phosphoric  acid  that 
has  low  vapor  pressure  is  used  as  an  electrolyte  to  suppress  the 
evaporation  and  dissipation  with  the  process  gases  (the  fuel  and 
oxidant  gases).  However,  the  decrease  in  cell  performance  as  a 
result  of  the  reduction  of  phosphoric  acid  is  not  negligible  at  the 
operating  temperature  of  about  200  °C.  For  stacks  that  are  designed 
to  have  an  operating  life  of  over  40,000  h,  a  basic  design  to  keep  low 
temperature  at  the  gas  exit  comparing  with  the  gas  inlet  of  the  cell 
is  adopted.  As  a  result  of  the  design,  the  evaporated  phosphoric  acid 
in  the  process  gas  is  condensed  and  recovered  to  the  cell  at  the  gas 
exit. 

Yoshioka  et  al.  have  conducted  experimental  and  numerical 
studies  for  evaporation  [1  ]  and  condensation  [2]  of  phosphoric  acid. 
In  those  studies,  the  evaporation  rate  was  treated  to  be  controlled 
by  the  diffusion  rate  in  the  concentration  boundary  layer  of  process 
gas  and  the  condensation  rate  was  treated  to  be  controlled  by  the 
cooling  rate  of  process  gas,  and  the  phosphoric  acid  distributions 
in  vapor  and  liquid  phases  and  the  recovery  of  phosphoric  acid  for 
the  model  cell  were  discussed. 

In  general,  the  vapor  pressure  of  phosphoric  acid  in  the  process 
gas  is  different  from  the  saturated  vapor  pressure,  as  the  evapora¬ 
tion  and  condensation  rates  between  the  cell  and  process  gas  have 
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limited  values.  Furthermore,  the  partial  pressure  of  phosphoric  acid 
in  the  process  gas  is  determined  by  complex  balance  of  evapora¬ 
tion  and  condensation  rates  and  passing  duration  time  of  process 
gas  that  go  through  the  cell  channel.  Therefore,  a  detailed  design 
method  that  includes  such  complex  mechanisms  was  required  to 
estimate  the  amount  of  the  dissipation  of  phosphoric  acid  precisely. 
Besides,  to  predict  the  phosphoric  acid  distribution  in  the  cell  with 
the  time  progress  is  important  to  estimate  the  degradation  of  cell 
performance. 

In  the  present  paper,  numerical  analysis  that  include  the  behav¬ 
ior  of  the  heat  transfer,  gas  flow,  and  evaporation  and  condensation 
of  phosphoric  acid  were  conducted  for  the  model  cell  that  con¬ 
sists  of  porous  electrode  substrate  and  matrix.  The  evaporation  and 
condensation  rate  of  phosphoric  acid  were  calculated  based  on  the 
nucleation  theory  for  a  two  components  system.  The  conditions  of 
the  temperature  and  impregnated  phosphoric  acid  concentration 
for  the  model  cell  were  correspondent  with  those  of  an  on  site  cell, 
and  time  series  solutions  for  the  duration  of 3300  h  were  calculated. 
The  phosphoric  acid  distributions  in  the  vapor  phase  and  liquid 
phase  at  the  process  gas  and  cell  were  calculated,  and  as  a  result 
of  those,  the  exhausting  and  remaining  phosphoric  acid  for  the  cell 
were  calculated,  and  the  results  were  evaluated  by  comparing  with 
the  results  for  the  model  cell  experiment. 

2.  Analysis  method 

The  analysis  code  that  was  developed  for  a  PAFC  stack  was  used 
for  the  calculations.  The  code  is  based  on  the  analysis  code  that 
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Nomenclature 

List  of  symbols 

J  evaporated  or  condensed  phosphoric  acid  flux 

(mols-1  cm-2) 

AG*  activation  energy  change  (J  mol-1  K-1 ) 
c  frequency  factor 

cE  frequency  factor  coefficient  for  evaporation 

cc  frequency  factor  coefficient  for  condensation 

Pa  partial  pressure  of  phosphorus  pentoxide 

Psa  saturated  vapor  pressure  of  phosphorus  pentoxide 
R  universal  gas  constant  (8.314 J  mol-1 I<-1 ) 

k  Boltzmann  constant  ( 1 .38E-23  J) 

t  temperature  (°C) 

T  absolute  temperature  (K) 

x  molar  fraction 

li  chemical  potential  (J  mol-1 ) 

y  activity  coefficient 

P  collision  rate 

Subscripts 

a  phosphorus  pentoxide 

w  water 

g  vapor  phase 

l  liquid  phase 

E  evaporation 

C  condensation 

s  saturation 

Superscript 

0  pure  component 


has  been  developed  for  the  characteristic  calculation  of  a  molten 
carbonate  fuel  cell  (MCFC)  stack  [3],  and  the  function  of  electrode 
reaction  has  modified  to  correspond  to  that  of  the  PAFC  and  the 
functions  for  the  evaporation  and  condensation  of  phosphoric  acid 
have  been  added.  Therefore,  the  analysis  features  about  the  heat 
transfer  and  gas  flow  are  the  same  as  those  of  the  code  for  MCFC. 

The  phosphoric  acid  migrations  in  PAFC  stack  are  considered  to 
be  caused  by  the  following  two  mechanisms. 

(1)  Migration  accompanied  with  the  phase  change  between  liq¬ 
uid  and  vapor.  That  is  dissipation  of  phosphoric  acid  from  the 
porous  electrode  to  process  gas  by  evaporation  and  recovery 
from  the  process  gas  to  porous  electrode  by  condensation  with 
the  transfer  by  gas  flow. 

(2)  Migration  in  the  liquid  phase.  That  is  migration  by  capillary 
force  caused  by  the  concentration  or  temperature  change  of 
phosphoric  acid  in  the  porous  electrode. 
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Fig.  1.  Outline  of  model  cell  experiment. 


surface,  R  is  the  universal  gas  constant  8.314  J  mol-1  K-1,  T  is  the 
liquid  phase  absolute  temperature  K,  AG*  is  the  activation  energy 
J  mol-1,  and  c  is  the  frequency  factor. 

The  calculation  methods  of  evaporation  and  condensation  rates 
will  be  mentioned  in  Appendix  A.  Whereas  saturated  vapor  pres¬ 
sure  of  P4O10  and  H20  were  derived  by  the  vapor  liquid  equilibrium 
for  a  two  components  system,  and  the  activation  energy  were 
derived  from  the  temperature  and  composition  for  the  vapor  phase 
and  liquid  phase  of  phosphoric  acid  in  the  process  gas  and  cell  based 
on  the  free  energy  change  with  phase  change.  The  concentrations  of 
P4O10  and  H20  in  the  evaporating  or  condensing  phosphoric  acid 
were  also  determined  from  the  free  energy  changes.  The  depen¬ 
dencies  of  the  temperature  and  concentration  of  P4O10  and  H20  on 
the  frequency  factor  were  derived  from  the  Hertz-Knudsen  equa¬ 
tion,  and  the  frequency  factor  coefficients  were  determined  so  as 
to  meet  the  analysis  results  with  the  experimental  results,  as  it  will 
be  mentioned  later. 

Whether  evaporation  or  condensation  arises  on  the  electrode 
surface  was  determined  by  comparing  P4O10  partial  pressure  Pa  in 
vapor  phase  at  the  process  gas  with  P4O10  saturated  vapor  pressure 
Psa  that  was  determined  by  the  composition  and  temperature  in 
liquid  phase  at  the  electrode  as  follows: 

Pa  <  Psa>  Evaporation  (2) 

Pa  >  Psa ;  Condensation  (3) 

The  code  for  MCFC  is  an  analysis  program  that  performs  rep¬ 
etition  calculation  and  acquires  a  steady  state  solution  about  the 
phenomena  that  are  the  results  of  coupling  for  the  electrochem¬ 
ical  reaction,  heat  transfer,  and  gas  flow.  However,  the  code  for 
PAFC  was  required  to  calculate  the  time  variation  of  phosphoric 
acid  dissipation  under  the  various  temperature,  vapor  phase  com¬ 
position,  and  liquid  phase  composition.  Therefore,  the  unsteady 
analysis  loop  that  includes  the  analysis  loop  of  the  code  for  MCFC 
was  programmed,  and  the  evaporation  and  condensation  rates  and 
the  amount  of  phosphoric  acid  in  the  vapor  phase  and  liquid  phase 
at  the  process  gas  and  electrode  were  calculated  in  this  unsteady 
loop.  The  calculation  for  the  balance  of  phosphoric  acid  in  unsteady 
analysis  loop  will  be  mentioned  in  Appendix  B. 


In  this  paper,  the  phenomenon  above  (1 )  that  has  main  effect  for 
exhausting  of  phosphoric  acid  was  considered  in  the  calculation. 
The  phosphoric  acid  was  regarded  as  a  two  components  system  of 
phosphorus  pentoxide  (P4O10)  and  water  (H20),  and  the  evapora¬ 
tion  and  condensation  rates  of  phosphoric  acid  between  the  cell  and 
process  gas  were  considered  to  have  limited  values,  and  these  were 
expressed  by  the  following  activation  process  based  on  the  nucle- 
ation  (droplet  condensation)  theory  of  a  two  components  system 
in  Refs.  [8-10]: 


J  =  c  exp 


(1) 


Here,  J  is  the  amount  of  evaporating  or  condensing  phospho¬ 
ric  acid  mol  s-1  cm-2  per  unit  area  and  unit  time  on  the  electrode 


3.  Analysis  conditions 

Analysis  was  conducted  on  the  conditions  that  correspond  to 
the  following  model  cell  experiment.  The  outline  of  the  exper¬ 
iment  is  shown  in  Fig.  1,  and  the  experimental  conditions  are 
shown  in  Table  1 .  The  experiment  was  carried  out  using  the  porous 
electrode  substrate  and  matrix  in  that  phosphoric  acid  was  impreg¬ 
nated,  and  the  steam  humidified  N2  gas  that  was  substitute  process 
gas  was  run  through  the  electrode  gas  channel.  The  temperature 
in  each  electrode  section  was  controlled  by  an  electrical  heater 
to  represent  the  cell  temperature  distribution  along  the  oxidant 
gas  flow  for  an  on  site  stack.  The  phosphoric  acid  migration  as 
liquid  phase  between  the  electrode  sections  was  prevented  by 
“Teflon”  boards  between  each  section.  The  amount  of  the  remain- 
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Table  1 

Experiment  and  analysis  conditions. 


Pressure 

1  atm 

Gas  flow  rate 

N2  (3.5 1  min-1 )  +  60 °C  humidified 

Operation  duration 

3300  h 

Initial  phosphoric  acid  amount 

11.4cm3  (85  wt%  H3PO4,  25°C)/section 

Temperature  Inlet  Section 

Outlet 

1 

2  3  4  5 

205  °C 

205  °C  198°C  178°C  163°C 

ing  phosphoric  acid  in  each  electrode  section  was  derived  from 
the  quantitative  analysis  results  by  ICP-OES  (Inductively  Coupled 
Plasma  Optical  Emission  Spectrometer)  after  3300  h  elapsed  state. 
Each  section  was  divided  into  nine  elements,  and  analysis  was  con¬ 
ducted  adopting  one  dimensional  model  in  the  direction  of  the 
gas  flow.  The  frequency  factor  coefficients  for  evaporation  and  a 
condensation  rate  were  determined  so  as  to  meet  the  remaining 
phosphoric  acid  of  analysis  results  with  that  of  the  experimental 
results.  At  this  time,  the  evaporation  rate  frequency  factor  coeffi¬ 
cient  Ce  in  Eq.  (A.2.13)  was  determined  according  to  the  amount 
of  remaining  phosphoric  acid  in  the  1st  section  of  electrode,  and 
the  condensation  rate  frequency  factor  coefficient  cc  in  Eq.  (A.2.1 7) 
was  determined  according  to  the  amount  of  remaining  phospho¬ 
ric  acid  in  the  5th  section  of  electrode.  The  same  frequency  factor 
coefficient  values  were  adopted  for  the  2nd,  3rd,  and  4th  electrode 
sections  and  the  frequency  factor  coefficients  used  for  the  analysis 
were  as  follows: 

Ce  =  -3.1E38  [m  s-1  kg1/2J1/,2mol_1/2]  (4) 

Cc  =  3. IE  -  81  [m  s_1  kg1/2J1/2mor1/2Pa_1]  (5) 

4.  Analysis  results 

4.1.  Evaporation  and  condensation  rates  of  phosphoric  acid 

The  evaporation  and  condensation  rate  distributions  of  phos¬ 
phoric  acid  along  the  gas  flow  direction  for  the  initial  state,  2000  h 
elapsed  state,  and  3300  h  elapsed  state  are  shown  in  Fig.  2.  Here, 
the  negative  rate  shows  evaporation  and  positive  rate  shows  con¬ 
densation.  The  evaporation  arises  at  sections  1  and  2  (0.0-0.4m,  is 
the  distance  from  the  gas  inlet).  The  condensation  arises  at  section 
5  (0.8-1. 0  m,  is  the  distance  from  the  gas  inlet)  where  the  elec¬ 
trode  temperature  fell  steeply  comparing  with  the  upstream.  The 
peak  magnitude  of  evaporation  and  condensation  rate  fell  accord¬ 
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Fig.  2.  Evaporation  and  condensation  rates  along  the  substrate  ( minus :  evaporation, 
plus:  condensation). 


Fig.  3.  P4O10  molar  fractions  of  evaporating  and  condensing  phosphoric  acid  along 
the  substrate. 


ing  to  the  progress  of  time,  and  the  peak  position  of  evaporation 
and  condensation  rate  gradually  moved  to  the  downstream  of  gas 
flow. 

The  evaporation  rate  had  gradual  distributions  in  wide  area, 
whereas  the  condensation  rate  had  steep  distributions  in  narrow 
area.  These  were  caused  by  the  following  mechanism.  The  activa¬ 
tion  energy  for  evaporation  rate  was  determined  by  thick  P4O10 
concentration  in  the  liquid  phase  at  electrode  and  response  to  the 
gas  flow  change  came  out  slowly,  though  activation  energy  for 
condensation  rate  was  determined  by  thin  P4O10  partial  pressure 
in  the  vapor  phase  at  process  gas  and  response  to  the  gas  flow 
change  came  out  rapidly,  as  shown  in  Eqs.  (A.2.1 1)  and  (A.2.1 6) 
respectively.  These  relations  were  also  shown  in  the  P4O10  par¬ 
tial  pressure  distribution  that  gradually  increased  at  the  inlet  and 
steeply  decreased  at  the  outlet  for  each  time  in  Fig.  6. 

P4O10  molar  fraction  distribution  of  evaporating  and  condensing 
phosphoric  acid  composition  along  the  gas  flow  direction  is  shown 
in  Fig.  3.  P4O10  molar  fraction  of  evaporating  composition  at  sec¬ 
tions  1  and  2  fell  according  to  the  time  progress,  but  there  was  few 
change  in  P4O10  molar  fraction  from  0.746  at  initial  state  to  0.744 
at  3300  h  elapsed  state.  In  condensing  composition  at  section  5, 
P4O10  molar  fraction  increased  according  to  the  time  progress,  but 
the  change  was  also  a  few  and  it  was  from  0.746  at  initial  state  to 
0.752  at  3300  h  elapsed  state.  P4O10  concentration  of  evaporating 
and  condensing  compositions  that  were  determined  by  Eq.  (A.2.1 0) 
were  thicker  comparing  with  the  phosphoric  acid  molar  fraction 
0.0922  (85  wt%  H3PO4)  that  was  impregnated  in  the  electrode  at 
the  initial  state. 


4.2.  Amount  of  phosphoric  acid  in  liquid  phase  at  the  electrode 

P4O10  and  H20  molar  concentration  distributions  in  liquid  phase 
at  electrode  along  the  gas  flow  direction  are  shown  in  Fig.  4.  The 
molar  concentration  of  P4Oi0  and  H20  decreased  gradually  at  sec¬ 
tions  1  and  2  by  evaporation,  and  increases  steeply  at  section  5  by 
condensation  with  the  progress  of  time.  Since  P4O10  concentrations 
of  evaporating  and  condensing  phosphoric  acid  were  thicker  than 
that  of  impregnated  phosphoric  acid  in  the  electrode,  the  varia¬ 
tion  of  P4O10  concentration  in  liquid  phase  at  electrode  was  larger 
than  the  variation  of  H20  concentration  at  sections  1,2,  and  5.  The 
peak  position  of  evaporation  and  condensation  rate  moved  toward 
the  gas  downstream,  and  the  evaporation  and  condensation  rate 
decreased  with  the  time  progress  in  Fig.  2.  These  were  correspon¬ 
dent  with  the  decrease  and  increase  of  the  molar  concentration  of 
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P4O10  and  H20  in  the  liquid  phase  according  to  the  gas  flow  for  the 
evaporation  and  condensation  respectively  in  Fig.  4. 

The  amount  of  remaining  P4O10  in  each  electrode  section  at 
3300  h  elapsed  state  for  the  analysis  results  and  experimental 
results  are  shown  in  Fig.  5.  As  mentioned  in  the  chapter  3,  the  fre¬ 
quency  factor  coefficient  of  evaporation  rate  was  determined  so  as 
to  meet  the  remaining  P4O10  of  analysis  result  in  section  1  with 
that  of  experimental  results,  and  the  frequency  factor  coefficient 
of  condensation  rate  was  determined  so  as  to  meet  the  remain¬ 
ing  P4O10  of  analysis  result  in  section  5  with  that  of  experimental 
results.  Regarding  the  frequency  factor  coefficients  for  sections  2-4, 
the  same  frequency  factor  coefficients  of  sections  1  and  5  were 
adopted,  and  as  a  result  of  those,  the  amount  of  remaining  P4O10  for 
sections  2-4  were  calculated  with  good  meet  to  the  experimental 
results.  Therefore,  adopting  appropriate  frequency  factor  coeffi¬ 
cients  for  the  evaporation  and  condensation  rates,  the  phosphoric 
acid  dissipation  for  other  temperature  and  gas  flow  conditions  was 
estimated  to  be  calculated  accurately  by  the  analysis  method  that 
mentioned  in  the  present  paper. 

The  amount  of  remaining  P4O10  for  the  analysis  and  experimen¬ 
tal  results  increased  according  to  the  section  goes  to  the  gas  outlet 
direction.  Especially  in  sections  4  and  5,  there  were  much  P4O10 
comparing  with  the  impregnated  P4Oi0,  and  the  recovery  of  phos¬ 
phoric  acid  from  the  process  gas  by  condensation  was  observed. 
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4.3.  Partial  pressure  and  saturated  vapor  pressure  of  phosphoric 
acid 

P4O10  partial  pressure  and  FI2O  partial  pressure  distributions  of 
vapor  phase  in  the  process  gas  along  the  gas  flow  direction  at  initial 
state,  2000  h  and  3300  h  elapsed  states  are  shown  in  Figs.  6  and  7 
respectively. 

At  the  initial  state,  the  partial  pressure  of  P4O10  and  H20 
increased  along  the  gas  flow  and  reached  to  saturation  in  section  1 
where  evaporation  occurred,  and  decreased  along  the  gas  flow  and 
reached  to  fixed  values  in  sections  4  and  5  where  condensation 
occurred.  The  gradual  increases  at  section  1  and  the  steep  decrease 
at  section  5  for  P4O10  and  H20  partial  pressures  were  correspon¬ 
dent  to  the  gradual  decrease  of  the  evaporation  rate  at  section  1 
and  the  steep  increase  of  condensation  rate  at  section  5  in  Fig.  2. 
The  increase  of  P4O10  and  H20  partial  pressures  at  the  gas  inlet  part 
became  loose  according  to  the  progress  of  time,  and  the  distance 
between  the  gas  inlet  and  the  position  where  the  vapor  pressure 
reached  to  saturation  became  long.  P4O10  partial  pressure  at  the 
gas  exit  that  determines  the  amount  of  exhausting  phosphoric  acid 
was  almost  stable  around  6.03E-3  Pa  in  the  progress  of  time.  H20 
partial  pressure  slightly  increased  from  6.88E-4  Pa  at  initial  state 
to  7.73E-4  Pa  at  3300  h  elapsed  state. 

P4O10  partial  pressure  in  the  vapor  phase  and  P4O10  saturated 
vapor  pressure  that  was  determined  by  the  condition  in  liquid 
phase  at  electrode  at  initial  state  and  3300  h  elapsed  state  are 
shown  in  Figs.  8  and  9  respectively. 


Fig.  5.  Residual  P4O10  in  the  electrode  after  3300  h. 


Fig.  7.  H20  partial  pressure  change  along  the  substrate. 
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Fig.  8.  P4O10  saturated  vapor  pressure  and  partial  pressure  along  the  substrate  at 
initial  state. 


At  the  initial  state,  P4O10  partial  pressure  at  the  gas  inlet  part 
was  lower  than  the  saturated  vapor  pressure  and  according  to  the 
gas  flow,  after  that  the  partial  pressure  increased  and  reached  to 
saturation.  At  sections  3  and  4  where  the  distance  from  the  gas 
inlet  part  are  from  0.4  to  0.8  m,  P4O10  partial  pressure  was  higher 
than  the  saturated  vapor  pressure  and  the  vapor  phase  of  process 
gas  was  supersaturating  state  regarding  to  the  temperature  and 
phosphoric  acid  concentration  in  the  liquid  phase  at  the  electrode. 
At  section  5  where  the  distance  from  the  gas  inlet  is  0.8  m  or  more, 
P4O10  partial  pressure  decreased  steeply  and  reached  to  saturated 
vapor  pressure  according  to  the  liquid  phase  conditions  at  the  gas 
exit. 

At  the  3300  h  elapsed  state,  P4O10  saturated  vapor  pressure 
increased  from  0  according  to  the  gas  flow,  and  P4O10  partial  pres¬ 
sure  distribution  was  almost  equal  to  the  saturated  vapor  pressure 
distribution  in  sections  1  and  2.  P4O10  saturated  vapor  pressure 
decreased  comparing  to  the  initial  state,  as  the  liquid  phase  phos¬ 
phoric  acid  concentration  was  decreased  because  of  the  thicker 
evaporating  P4O10  concentration  as  shown  in  Fig.  3.  P4O10  par¬ 
tial  pressure  was  supersaturating  state  at  sections  3  and  4,  and  it 
decreased  steeply  at  section  5  and  reached  to  the  saturated  vapor 
pressure  at  the  gas  exit.  The  saturated  vapor  pressure  steeply  rose 
at  the  gas  inlet  portion  of  section  5  caused  by  the  steep  increase  of 
P4O10  concentration  in  the  electrode  by  condensation  as  shown  in 
Fig.  4. 


0  0.2  0.4  0.6  0.8  1 

Distance  from  Inlet  (m) 

Fig.  9.  P4O10  saturated  vapor  pressure  and  partial  pressure  along  the  substrate  after 
3300  h. 
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Fig.  10.  Time  variation  of  evaporation  and  condensation  rates  (minus:  evaporation, 
plus:  condensation). 
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Fig.  11.  Time  variation  of  P4O10  partial  pressure. 


4.4.  Variation  in  the  time  progress 

The  phosphoric  acid  evaporation  and  condensation  rates  at  each 
section,  P4O10  partial  pressures  in  the  vapor  phase  at  process  gas, 
the  amounts  of  remaining  P4O10  in  the  liquid  phase  at  electrode 
at  each  section  for  the  variation  with  time  progress  are  shown  in 
Figs.  10-12  respectively. 

Regarding  the  phosphoric  acid  evaporation  and  condensation 
rates,  the  average  evaporation  rate  at  section  1  decreased  and 
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Fig.  12.  Time  variation  of  residual  P4O10  in  the  electrode. 
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the  average  evaporation  rate  at  section  2  increased  with  the 
progress  of  time.  These  were  caused  by  the  difference  between 
the  saturated  vapor  pressure  and  the  partial  pressure  of  P4O10 
were  decreased  at  section  1  and  increased  at  section  2  with  the 
progress  of  time  as  shown  in  Figs.  8  and  9  respectively.  The 
evaporation  and  condensation  rate  at  section  3  was  maintained 
almost  the  same  as  initial  state,  and  the  average  condensation 
rate  at  sections  4  and  5  decreased  slightly  with  the  progress 
of  time. 

P4O10  partial  pressure  at  sections  1  and  2  decreased  with  the 
progress  of  time  according  to  the  decrease  of  P4O10  amount  in  the 
liquid  phase  at  electrode  as  shown  in  Fig.  12.  P4O10  partial  pressure 
at  sections  3  and  4  was  maintained  almost  the  same  until  2000  h 
and  after  that  decreased  slightly,  these  were  caused  by  P4O10  sat¬ 
urated  vapor  pressure  decrease  at  3300  h  elapsed  state  comparing 
with  the  initial  state  as  shown  in  Figs.  8  and  9.  In  section  5,  P4O10 
partial  pressure  increased  with  the  time  progress  caused  by  the 
increase  of  P4O10  amount  in  the  liquid  phase  at  electrode  as  shown 
in  Fig.  4. 

The  amount  of  P4O10  reduced  0.0322  mol  and  0.00930  mol  in 
sections  1  and  2  respectively  at  the  3300  h  elapsed  state  compar¬ 
ing  with  the  initial  state.  P4O10  at  section  3  was  maintained  around 
0.041 9  mol  through  the  duration  of  3300  h.  The  amount  of  P4O10 
increased  0.0125  mol  at  section  4  and  0.0274  mol  at  section  5  at 
the  3300  h  elapsed  state  comparing  with  the  initial  state.  Conse¬ 
quently,  the  P4O10  at  sections  1  and  2  migrated  to  sections  4  and  5 
by  the  evaporation  and  condensation  accompanied  with  gas  flow. 
The  total  remaining  P4O10  in  the  electrode  reduced  0.00210  mol 
after  3300  h,  this  shows  the  amount  of  exhausted  P4O10  with  the 
process  gas  and  it  was  1.0%  of  the  total  impregnated  P4O10  at  the 
initial  state. 


5.  Conclusion 

The  evaporative  dissipation  and  migration  characteristics  of 
phosphoric  acid  in  the  model  cell  were  calculated  by  the  analy¬ 
sis  code  for  PAFC  that  include  the  behavior  of  the  heat  transfer, 
process  gas  flow,  and  evaporation  and  condensation  of  phos¬ 
phoric  acid.  The  analysis  was  conducted  for  the  model  cell  that 
consists  of  porous  electrode  substrate  and  matrix  and  has  corre¬ 
spondent  conditions  with  that  for  an  on  site  stack  in  the  duration 
of  3300  h. 

The  evaporation  and  condensation  behavior  of  phosphoric  acid 
were  calculated  based  on  the  nucleation  theory  of  a  two  compo¬ 
nents  system,  and  consequently  the  phosphoric  acid  dissipation 
characteristic  for  each  electrode  section  was  calculated.  The  analy¬ 
sis  results  had  good  meet  with  the  experimental  results  by  adopting 
appropriate  frequency  factor  coefficients  for  the  evaporation  and 
condensation  rates.  Consequently  the  effectiveness  of  the  present 
analysis  method  was  validated. 

P4O10  partial  pressure  was  affected  by  the  conditions  of  the 
vapor  phase  in  process  gas  and  the  liquid  phase  in  electrode,  and 
for  some  conditions,  P4O10  partial  pressure  could  be  higher  than 
the  saturated  vapor  pressure  according  to  the  liquid  phase  in  the 
electrode,  and  the  supersaturating  states  were  resulted  to  be  pos¬ 
sible. 

The  amount  of  remaining  P4O10  decreased  in  the  high  tempera¬ 
ture  sections  at  gas  inlet  side  and  increased  in  the  low  temperature 
sections  at  gas  outlet  side  through  the  3300  h  duration.  Con¬ 
sequently,  P4O10  in  the  gas  inlet  sections  migrated  to  the  gas 
outlet  sections  by  the  evaporation  and  condensation  accompa¬ 
nied  with  the  gas  flow.  Besides,  1.0%  of  total  impregnated  P4O10 
was  exhausted  with  process  gas  after  3300  h.  The  remaining  dis¬ 
tribution  and  exhausting  amount  of  P4Oi0  for  the  electrode  were 
estimated  into  the  details. 


Appendix  A.  A.  Evaporation  and  condensation  rates  of 
phosphoric  acid 

AA.  The  vapor  liquid  equilibrium  of  a  two  components  system 


The  phosphoric  acid  was  regarded  as  a  two  components  system 
of  phosphorus  pentoxide  (P4O10)  and  water  (FI20),  and  the  vapor 
liquid  equilibrium  was  considered.  The  saturated  vapor  pressure 
Psa  for  P4O10  in  arbitrary  concentration  of  phosphoric  acid  in  the 
liquid  phase  was  derived  as  follows: 


Psa  =  P&-Ya-Xa  (A.1.1) 

Here,  P°Q  is  the  pure  P4O10  saturated  vapor  pressure  mmHg, 
ya  and  xa  are  the  activity  coefficient  and  the  molar  fractions  of 
P4O10  in  the  phosphoric  acid  respectively.  P°a  was  derived  from 
the  Antoine  equation  [4],  and  ya  was  derived  by  regressing  P4O10 
vapor  concentration  data  shown  in  [5]  to  the  equation  shown  in  [6] 
as  follows: 

logic  Ps°a=A— ^  (A.1.2) 


Ya  =  exp 


fla(l  —  *a)  +  ba(  1 
RTi 


(A.1.3) 


Flere,  t\  is  the  liquid  temperature  °C,  R  is  the  universal  gas  con¬ 
stant  8.314 J  mol-1  K-1,  T/  is  the  liquid  absolute  temperature  (I<), 
and  each  coefficient  was  as  follows: 


A  =  9.7070,  B  =  3822,  C  =  201.0 


aa  =  — 6.732E4,  ba  =  -4.607E4 


Flere,  aa  and  ba  were  regressed  by  the  almost  same  coefficient 
in  the  temperature  range  of  Tt  =  436-491  K. 

Similarly,  the  saturated  vapor  pressure  Psw  for  H20  in  arbitrary 
concentration  of  phosphoric  acid  in  the  liquid  phase  is  derived  as 
follows: 


Psw  —  Psw  '  Yw  •  (1  —  *a)  (A.1.4) 

Flere,  P°w  is  the  pure  H20  saturated  vapor  pressure  mmFIg,  yw  is 
the  activity  coefficient  of  H20,  and  xa  is  the  molar  fractions  of  P4O10 
in  phosphoric  acid.  Ps°w  was  derived  from  the  Wagner  equation  [7], 
and  yw  was  derived  by  regressing  the  saturated  vapor  pressure  data 
of  FI20  in  phosphoric  acid  shown  in  [5]  to  the  equation  shown  in 
[6]  as  follows: 

In  =(a&  +  bV'-5  +  c&3+d#6)  jT^,  &  =  \-T 

C  (A.1.5) 


Yw  =  exp 


awXa  +  bwXa  \ 

PT[  j 


(A.1.6) 


Flere,  the  critical  pressure  Pc  was  1 .6591 E5  mmFIg,  and  the  criti¬ 
cal  temperature  Tc  was  647.3  K  and  each  coefficient  was  as  follows: 


a  =  -7.76451,  b  =  1.45838,  c  =  -2.77580,  d  =  -1.23303 


Coefficients  aw  and  bw  vary  with  temperature,  and  were 
regressed  as  follows  about  each  temperature: 


Ti  =  373  K  : 
Tt  =  398  K  : 
Tt  =  423  K  : 
1/  =  448  K  : 


aw  =  5.499E5,  bw  =  -1.278E6 
aw  =  — 4.904E5,  bw  =  -1 .695E6 


aw  =  5.538E5,  bw  =  -  1.051E6 


aw  =  -  5.229E4,  bw  =  -5.117E6 
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T,  =  473  K 
Ti  =  498  K 
1/  =  523  K 

A.2.  Evaporation  and  condensation  rates 


aw  =  — 5.454E4,  bw  =  -5.263E6 
aw  =  -1.230E6,  bw  =  2.660E6 
aw  =  — 7.305E5,  bw  =  -3.496E5 


The  evaporation  and  condensation  rates  of  phosphoric  acid  were 
expressed  as  following  activation  process,  based  on  the  theory  of 
the  nucleation  (droplet  condensation)  of  a  two  components  system 
mentioned  in  Refs.  [8-10]: 


J  =  c  exp 


-AG* 

~rT 


(A.2.1) 


Here  J  is  the  amount  of  phosphoric  acid  that  evaporate  or  con¬ 
densate  per  the  unit  area  on  an  electrode  surface  and  unit  time 
mol  s-1  cm-2,  R  is  the  universal  gas  constant  8.314J  mol-1 1<-1 ,  and 
J  was  treated  as  minus  and  plus  for  evaporation  and  condensation 
respectively,  T  is  the  liquid  absolute  temperature  (K),  AG*  is  the 
activation  energy  (J  mol-1 ),  c  is  the  frequency  factor. 

A.2.1.  Evaporation  rate 

The  free  energy  change  A GEa  accompanied  with  the  phase 
change  in  evaporation  was  derived  as  follows: 

A Gp  ( Xa ,  Xw)  =  Xa A/Xa  +  Xw A flw  =  Xaiflai  —  P'ag)  +  Xw(Mw/  —  P'Wg) 

(A.2.2) 


=  xaRT  In  ^.+xwRT  ln^A 

1  sa  1  sw 


Here,  xa  and  xw  are  the  molar  fraction  of  P4O10  and  H20  in  evap¬ 
orating  liquid  phase,  A fia  and  A/xw  are  chemical  potential  change 
(Jmol-1)  accompanied  with  the  phase  change  of  P4Oi0  and  H20. 
Psa  and  Psw  are  the  saturated  vapor  pressure  of  P4O10  and  H20  in 
the  liquid  phase  at  electrode  respectively,  Psa  and  Psw  are  saturated 
vapor  pressure  of  P4O10  and  H20  for  the  liquid  phase  correspond  to 
evaporating  vapor.  T is  liquid  absolute  temperature  at  electrode  (I<). 

The  following  relation  is  satisfied  in  the  liquid  phase  composi¬ 
tions  for  evaporating  vapor: 

xw  =  l— xa  (A.2.3) 

The  following  relation  was  derived  from  Eqs.  (A.2.2)  and  (A.2.3) 


A  CE(xa)  =  xaRT  In  ^  +  (1  -  xa)Rrin  ^ 

1  sa  1  sw 


(A.2.4) 

Psa*  Psw  and  Psa*  Psw  were  derived  from  Eqs.  (A.1.1)  and  (A.1.4) 


Psa(Xa,  T )  =  Psa(T)  '  *a  •  Ya(Xa,  T) 

(A.2.5) 

Psw(Xw,  T)  =  P$w(T )  '  •  Yw{Xw,  T) 

(A.2.6) 

Psa(Xa,  T)  =  Psa(P)  '  Xa  •  Yai^a,  T) 

(A.2.7) 

Psw{Xw,  T )  =  PS°W(T)  •  Xw  •  Yw(Xw,  T) 

(A.2.8) 

Here,  xa  and  xw  are  molar  fraction  of  P4O10  and  H20  of  liquid  at 
electrode. 

Eq.  (A.2.4)  is  transformed  as  follows  using  the  equations  from 
(A.2.5)  to  (A.2.8),  (A.1.3)  and  (A.1.6): 

A CE(xa)  =  XaRT  In  ^4^1  S  +  (1  -  Xa)RT  In  -  *wKw(*w’ 


=  xaRT  In 


xaya(xa,  T ) 

XaYa{Xa ,  T) 


(1  —  Xa)Yw(l  —  Xa,  T ) 


xa  exp((a„(l  -Xa)2  +ba(  1  -Xa)3)/RT) 


+  (1  -Xa)RT  In 


xwYw(xw,  T ) 


(1  -xa)exp((awx2  +bwxl)/RT) 


(A.2.9) 


In  Refs.  [8-10],  the  phosphoric  acid  composition  of  condensa¬ 
tion  was  derived  as  a  result  of  searching  a  saddle  point  for  the  free 
energy  change  that  is  equivalent  to  Eq.  (A.2.2).  The  evaporation 
composition  was  derived  at  the  critical  point  of  Eq.  (A.2.9)  based 
on  above. 

In  order  to  search  a  critical  point,  Eq.  (A.2.9)  was  differentiated 
by  xa  and  substituting  the  result  as  0,  the  following  equation  was 
derived: 


4  (ba  +  bw)Xa  +  3(— aa  —  3  ba  +  aw  —  bw)x 2  4-  2(2  aa 

+  3 ba  -  aw)xa  -  aa  -  ba  +  RT  In  —  ~-*a)*aKa  =0 

xaxw  Yw 


(A.2.10) 


P4O10  molar  fraction  in  the  evaporated  phosphoric  acid  was 
derived  by  solving  the  above  equation  with  xa  numerically.  The 
activation  energy  A G|  was  derived  as  the  free  energy  change  on 
critical  point  by  substituting  x*  for  xa  in  Eq.  (A.2.9)  as  follows: 


A  G*e  =  AGf(x*)  =  x*RT  In 


+  (1  -  x*)RT  In 


XaYa{xa,  T ) 


x*exp((afl(l  -x*)  +jba(l  -X*)  )/RT) 


XwYw(xw,  T ) 


(1  -  x^)exp  ((dwXS  +  bwxf )  /RT) 


(A.2.1 1) 


Here,  the  free  energy  change  direction  was  determined  so  as  to 
satisfy  that  the  evaporation  direction  becomes  minus  for  J  in  Eq. 
(A.2.1),  and  as  a  result  of  this,  the  following  relation  was  satisfied: 

-AG*  =  AG*e  (A.2.12) 

The  frequency  factor  Ce  was  derived  as  follows: 

Ce  =  Ce  ■  PaE  ■  S  ■  cwt  (A.2.1 3) 

Here,  cE  is  frequency  factor  coefficient,  S  is  area  ratio  that  derived 
by  dividing  surface  area  of  electrode  by  appearance  surface  and 
treated  as  1.  cw/  is  H20  molar  concentration  in  liquid  phase  at  the 
electrode  mol  cm-3.  paE  is  collision  rate  of  P4Oi0  in  liquid  phase 
at  the  electrode  and  derived  assuming  Herz-Knudsen  equation  for 
condensation  [11]  as  follows: 


PaE  = 


27TMakT 


(A.2.14) 


Here,  xa  is  the  molar  fraction  of  P4O10  in  the  liquid  phase,  Ma 
is  the  molecular  weight  of  P4O10,  k  is  the  Boltzmann  constant 
1.38E— 23JK-1,  and  T  is  the  liquid  phase  absolute  temperature  (I<). 

Consequently,  the  evaporation  rate  is  derived  adopting  activa¬ 
tion  energy  (A.2.12)  and  frequency  factor  (A.2.13)  as  follows: 


/AG* 

JE  =  C£exp(  — 


(A.2.15) 


A.2.2.  Condensation  rate 

According  to  the  same  manner  for  the  evaporation  rate,  the  acti¬ 
vation  energy  AG* ,  frequency  factor  cc,  collision  rate  in  the  vapor 
phase  Pac *  and  condensation  rate  Jc  were  derived  as  follows: 


AG£  =  AGc  (x*)  =  -x*aRT  In 


-(1  -x*JRT  In 


psaxa exp((aa(l  -x*a)2  +  ba(  1  -x*f)/RT) 


0 1  -x*)exp  ((awxf  +bwxf)  / RT ) 


cc  =  CC  ■  PaC  •  5  •  Cwg 

a  _  Pa 

PaC  — 


^/2jrMJcT 


(A.2.16) 

(A.2.17) 

(A.2.18) 
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Jc  =  cc  exp 


(A.2.19) 


Here,  Pa  and  Pw  are  the  partial  pressure  of  P4O10  and  H20  in 
condensing  vapor  phase  phosphoric  acid;  cc  is  the  frequency  factor 
coefficient,  cwg  is  the  molar  concentration  of  H20  in  vapor  phase 
mol  cm-3,  paC  is  the  collision  rate  of  P4O10  in  the  vapor  phase  at 
process  gas. 


Appendix  B.  B.  Unsteady  phosphoric  acid  balance 

The  balance  of  phosphoric  acid  between  the  liquid  phase  in  elec¬ 
trode  and  the  vapor  phase  in  process  gases  were  regarded  as  the 
P4O10  balance  and  H20  balance  separately. 

The  balance  for  P4O10  was  derived  as  follows: 

v^f=AJx*a  (B.l) 

Qin  •  cagin  =  A  J  •  Xq  +  qout  ■  cagout  (B.2) 

Here,  V  and  A  are  the  volume  (cm3)  and  surface  area  (cm2)  of 
the  electrode  element  respectively.  J  is  the  evaporation  or  con¬ 
densation  rate  of  phosphoric  acid  (mol  cm-2  s-1),  x*  is  the  molar 
fraction  of  P4O10  in  the  evaporating  or  condensing  phospho¬ 
ric  acid  vapor,  cai  is  the  molar  concentration  of  P4O10  in  liquid 
phase  phosphoric  acid  at  electrode.  The  qin  and  qout  are  the  vol¬ 
ume  flow  rate  of  the  process  gases  at  the  inlet  and  outlet  of  the 
electrode  element  channel  respectively  (cm3s-1),  cagin  and  cagout 
are  the  molar  concentration  of  P4O10  in  vapor  phase  phosphoric 
acid  at  the  inlet  and  outlet  on  the  element  channel  respectively 
(mol  cm-3). 


The  balance  for  H20  is  derived  as  following: 

V^L=A-J-(  1-xS)  (B.3) 

Qin  •  cwgin  =  A  J  ■  ( 1  —  Xq)  +  q0ut  •  Cwgout  (B.4) 

Here,  cw/  is  H20  molar  concentration  (mol  cm-3 )  in  liquid  phase 
phosphoric  acid  at  electrode;  cwgin  and  cwgout  are  the  molar  concen¬ 
tration  of  H20  in  the  vapor  phase  phosphoric  acid  at  the  inlet  and 
outlet  of  the  element  channel  respectively  (mol  cm-3). 

The  difference  equations  for  the  equations  from  (B.l  )-(B.4)  were 
calculated  at  the  unsteady  loop  that  including  the  electrochem¬ 
ical  reaction,  heat  transfer,  and  gas  flow  also,  and  consequently 
the  molar  amount  of  P4O10  and  H20  in  the  electrode  element  in 
each  time  step  were  calculated  as  Vx  cai  and  V x  cwi  respectively. 
The  mol  concentration  of  P4O10  and  H20  in  the  vapor  phase  pro¬ 
cess  gases  in  each  time  step  were  calculated  as  ( cagin  +  cagout)l 2  and 
( cwgin  +  cwgout)l 2  respectively. 
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1  Product  names  mentioned  herein  may  be  trademarks  of  their  respective  com¬ 
panies. 


